Shot noise of interference between independent atomic systems. 
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We study shot (counting) noise of the amplitude of interference between independent atomic 
systems. In particular, for the two interfering systems the variance of the fringe amplitude decreases 
as the inverse power of the number of particles per system with the coefficient being a non-universal 
number. This number depends on the details of the initial state of each system so that the shot noise 
measurements can be used to distinguish between such states. We explicitly evaluate this coefficient 
for the two cases of the interference between bosons in number states and in broken symmetry states. 
We generalize our analysis to the interference of multiple independent atomic systems. We show 
that the variance of the interference contrast vanishes as the inverse power of the number of the 
interfering systems. This result, implying high signal to noise ratio in the interference experiments, 
holds both for bosons and for fermions. 



Interference between independent systems is an in- 
teresting phenomenon, which goes back to the discov- 
ery of the Hanbury-Brown-Twiss (HBT) effect (see e.g. 
Ref. The very possibility that particles from inde- 
pendent condensates can interfere, i.e. have a certain 
relative phase, is quite intriguing and is not entirely ob- 
vious Q . The origin of this phenomenon is the quantum 
indistinguishability of identical particles and the super- 
position principle. Recent experimental advances in cold 
atom systems opened new possibilities to detect the inter- 
ference both in equilibrium 0, [3, H, EJ and in nonequilib- 
rium [?[ situations. Such experiments even allow probing 
various properties of interacting systems. For example, 
an analysis of the scaling of the interference signal with 
the imaging size lead to the direct observation of the 
Kosterlitz-Thouless phase transition Q . The problem of 
interference also attracted a lot of theoretical attention. 
In Ref. Q Javanainen and Yoo numerically studied the 
interference between two Fock states for a given set of 
detectors. It was shown that the resulting count dis- 
tribution of atoms was similar to the one arising from 
the interference of two condensates with randomly bro- 
ken phases. Later Castin and Dalibard analyzing another 
Gedanken experiment came to the same conclusion and 
showed that these random phases are induced by the de- 
tectors 0]. Emerging interference between two number 
states was later interpreted as an indication of the pres- 
ence of hidden variables [xoj | . In Refs. 11, 12 1 the authors 
studied the interference between independent fluctuating 
condensates and reached a similar conclusion that for the 
large ideal condensates the amplitude of interference does 
not fluctuate. Recently Altman et. al. suggested using 
noise interferometry as a new probe of interacting atomic 
systems [13 1. These ideas were later experimentally im- 
plemented 13, Ell- 
in fact the reason why two independent condensates 
with fixed number of particles in each must interfere fol- 
lows from the basic principles of quantum mechanics. 
Indeed, because of the number-phase uncertainty each 
of the condensates does not have a well defined phase. 
However, according to quantum mechanics any measure- 



ment probe sensitive to the phase difference (whether it 
is a time of flight image or something else) will project 
the condensates to the state with a well defined relative 
phase. As a result the condensates will interfere but the 
relative phase will be random for each experimental run 
in agreement with Refs. EH as well as this work. 

The main purpose of this paper is to study shot noise 
of the interference between independent atomic systems. 
For the most part we will consider bosonic systems (con- 
densates for short). However, our approach is indepen- 
dent of the atom statistics and where necessary we will 
give the explicit expressions for fermions as well. 

We point that if the systems are independent then 
the interference between them can be considered as a 
noise, similarly to the HBT effect. Indeed, the average of 
the atom density over many experimental runs does not 
result in any oscillating component. However, for the 
two condensates with large number of particles, each run 
yields well defined interference fringes 
fringe amplitude is thus a quantity which does not aver- 
age to zero. Moreover, as we mentioned above, this am- 
plitude is expected to have vanishing fluctuations. Con- 
versely, if the number of particles in each of the two sys- 
tems is small, we expect to see significant shot noise and 
large fluctuations of the fringe amplitude. The analysis 
of these fluctuations is the subject of the present inves- 
tigation. On passing on, we mention that shot noise ex- 
periments are a very powerful tool in condensed matter 
physics (see Ref. [16j for a review) and in quantum optics 
(see e.g. Ref. El). Such experiments can be used to de- 
tect charge of quasi-particles, transmission properties of 
small conductors, entanglement between electrons, allow 
to distinguish quantum and classical light sources, etc. 

We emphasize that when one analyzes the interference 
between extended systems, one encounters at least two 
different sources of noise (apart from the noise associated 
with the probing beam). One of them originates from the 
phase fluctuations within each system, which can have 
either quantum or thermal ori gin . This type of noise 
was analyzed in detail in Refs. [Ill E^. In particular, it 
was shown that the scaling of the fringe amplitude with 
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the system size and its distribution function contain in- 
formation about phase correlation functions within each 
system. The second source of noise, which is studied in 
the present paper, has a purely quantum nature origi- 
nating from the commutation relations between identical 
particles. This counting or shot noise exists even in ideal 
nonintcracting systems and as we show below it contains 
important information about the nature of the interfering 
states. 
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FIG. 1: Schematic view of the interference experiment be- 
tween two independent systems. 



Assume that all bosons in each of the two interfer- 
ing systems systems (1 and 2) occupy the same quantum 
state. We consider the standard time of flight measure- 
ment setup (see Fig. [IJ. At a certain moment in time 
these particles are allowed to expand. For simplicity we 
assume that the particles expand only in one direction x. 
For the two dimensional systems with strong transverse 
confinement like in e.g. Ref. @ this assumption is usu- 
ally well justified since the transverse kinetic energy is 
very large. For the system of two one-dimensional con- 
densates 0| the expansion for each condensate occurs in 
the radial direction . In this case one can always ana- 
lyze the image in the plane of the two condensates and 
the present analysis remains intact. In a more general 
case one has to compute the overlap between the expand- 
ing Wannier functions corresponding to the two different 
condensates. However, we emphasize that this overlap 
brings a more complicated geometrical factor character- 
izing the interference, which enters only as the prefactor 
into the interference amplitude and does not affect main 
conclusions of this paper. After waiting until the size 
of the clouds greatly exceeds the original separation be- 
tween the condensates d, the density of these atoms is 
measured using probing beams. We assume that there is 
no photon shot noise and the atoms are detected with a 
100% probability. For simplicity we also assume that d is 
much larger than the initial size of each condensate wq. 

As we already argued, there is no average interference 
contrast between the two condensates. However, there is 
a well defined observable A 2 , equal to the square of the 
interference amplitude [ill ]. Let us define the fluctuating 



variable A, average of which over many experimental runs 
gives A2, in the following way: 



A = 



dxdx'p(x, t)p(x', t)e lQ 



(x— X ) 



dxp(x, t), 



(1) 

Here p{x 1 1) is the number of the absorbed photons at 
the time t at the position x and Q = md/Ht, where m is 
the atom's mass. The integral over x and x' should be 
understood as the sum, where the discretization step is 
determined by the detector and can not be smaller than 
the photon wavelength. Apart from the second term in 
Eq. ([1]), which as we will see below has a quantum ori- 
gin, the expression for A is just a square of the Fourier 
transform of the absorbtion image taken at the wave vec- 
tor Q. Note that since by assumption all the detectors 
have 100% efficiency, the number of the absorbed pho- 
tons coincides with the number of atoms at the detection 
point: p(x,t) = n(x,t) = aJ (x,t)a(x,t), where a'(x,t) 
and a(x, t) are the time-dependent creation and annihi- 
lation operators of the atoms. We understand the equal- 
ity sign here and in expressions below in the sense that 
the statistical properties of p{x, t) and the corresponding 
quantum operator are equivalent. 

Using the bosonic commutation relations 
[a(x, t),a^{x',t)\ = 8{x — x') it is easy to see that 
Eq. |T]) can be rewritten as 
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dxdx'aJ (x, t)a'(x' , t)a(x, t)a(x' ,t)e 



(2) 

We can further simplify Eq. ([2]) using that in the long 
time limit 0, E3] 

a(x,t) ^ a lU (x,t)e iQl ^ x + a 2 u{x,t)e iQ2{t)x , (3) 

where Q12 = mix ± d/2)/ht, 0,1.2 are the bosonic oper- 
ators in the Schrodinger representation corresponding to 
the systems 1 and 2, and u(x,t) is the time-dependent 
Wannier function corresponding to either of the two con- 
densates (for simplicity we assume that these Wannier 
functions are identical). Using Eqs. (J2J> and ([3|) we find 



A 9 =A 



{a\a\a2(ii) 



(4) 



where the overline implies averaging over many experi- 
mental runs and the angular brackets denote the expec- 
tation value. We note that if the initial width of the 
condensates wq is finite then there will be corrections to 
Eq. (|4|) proportional to exp(— d 2 /wq). However, they are 
negligible as long as d ^> wo. 

We would like to emphasize that it is crucial to first 
take the normal order of the operator A (or any other 
operator), like we did in Eq. ([2]), and only then use 
Eq. ([3|). Indeed, it is easy to check that the asymptotic 
form ([3|) does not preserve the correct commutation re- 
lations between the operators a(x,t) and a)(x' ,t). It is 
therefore very important to order the operators first and 
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use Eqs. ([4} only after that. It is straightforward to check 
that if one proceeds in a different way then one will get 
un-physical results with negative probabilities for the ob- 
servables. 

One can similarly analyze the fermionic case. The re- 
sulting expression for the interference amplitude is iden- 
tical to Eq. (|4]). There is, however, one subtlety. Namely 
the mean value of A defined in Eq. ([1]) is negative, reflect- 
ing the anti- bunching of fermions. Therefore the ampli- 



tude of interference in the fermionic case should be de- 
fined as A 2 — — A. 

Next we want to find the fluctuations of A. For this 
purpose we need to compute A4 squaring Eq. ([TJ) and 
substituting p(x,t) by n(x,t). As we mentioned above 
we need to normal order the resulting expression first and 
then use the asymptotic forms of the operators a{x, t) and 
a){x,t). As a result one finds 



A 4 = J J J J dx 1 dx 2 dx' 1 dx' 2 (: n(x 1 ,t)n(x 2 ,t)n{x' 1 ,t)n(x , 2 ,t) : ) e *Q(*i+**-*'i-*ai 

+ 2JJ Jdxidx 2 dx' 1 (:n(x 1 ,t)n(x 2 ,t)n{x' 1 ,t):)e iQ( - X2 ~ x i ) +J J dx x dx % (; n(x 1 ,t)n(x 2 ,t) :) , (5) 

where the semicolons imply the normal-ordered form, i.e. that the creation operators appear on the left of the 
annihilation operators. Substituting Q into the integrals above yields 



A4 = (a\ 2 a\ 2 a 2 a 2 } + 2 (a[ 2 o4 a 2 a\) + 2 {a\a\ 2 a\ a\) + (a\ 2 a 2 ) + (a^af) + 2{a\a\a 2 ai) . 

I 



(6) 



Let us define the relative width of the distribution 
w = \J A 4 — A\ /A 2 . We explicitly look into two dif- 
ferent initial states. First we consider the interference 
between two independent coherent states with on aver- 
age N atoms in each state. In this case from Eqs. (0| 
and (HO) we find: 



A 2 = N 2 , A 4 = A 4 + 4A 3 + 4A 2 



2y/N + T 

N 



(7) 
(8) 



As expected, w vanishes as N — > 00. 

Next we consider the interference between the two 
number states with N atoms in each of them. Then 



A 2 = N 2 , A 4 = N 4 + 2N 3 + N 2 - 2N => 
V2 

w = —= 



1 



1 

2N 



1 

N*' 



(9) 
(10) 



Asymptotically w also decreases as the inverse square 
root of the number of particles at large N. However, the 
coefficient appears to be smaller by a factor of \p2 than 
in the case of the two coherent states. 

We note that Eq. © is also valid for fermions. How- 
ever, in this case all the terms except the last one iden- 
tically vanish and thus A4 — 2A\ independent of the 
details of the fermionic state. This result comes from 
the fact that in this simple setup there are at most two 
interfering particles. 

In the case of the interference of extended conden- 
sates, which exp and only in the transverse directions (see 
Refs. 0, [llL E3) ' t ne expressions above are easily gener- 
alized. For example in Eq. © one has to do the following 
substitutions: 



/ f2 t2 2 2\ 



I t 2 t 2\ 



(a\ 2 a 2 ) 



{a\a\a 2 a%) 



dzidz 2 dz 3 dz 4 {a\{zi)a\{z 2 )a\{zz)a\{z 4 )a 2 {zi)a 2 {z 2 )ai{z 3 )ai{z 4 )) , (11) 

dz x dz 2 dz 3 (a{(z 1 )a\(z 2 )al(z 3 )a 2 (z 2 )a 1 (z 3 )a 1 (z 1 )), (12) 

dzidz 2 (a\(zi)a\(z 2 )ai(z 2 )ai(zi)), (13) 

dz x dz 2 (a\(z 1 )al(z 2 )a 2 (z 1 )a 1 (z 2 )) . (14) 



Other substitutions can be obtained from these expres- sions by simple permutations. In the equations above Zj 
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denote coordinates along the condensates and the inte- 
gration is taken over the area or the length of the con- 
densates. For the systems with long-range correlations 
the terms containing the largest number of bosonic op- 
erators give the leading contribution to A4 in agreement 
with the results of Refs. 11, [l^. Indeed, for example, 



for the interference between two ID condensates at zero 
temperature Eqs. ((U and (fTT|) - (fl4|) give 

-y K (l + B(pL) 1 / K - 1 ), (15) 



A A = AL 



where A and B are the non-universal constants, p is the 
mean particle density, L is the system size or the imaging 
length (see Refs. [HHllI for more details). The Luttinger 
parameter K above characterizes the strength of the in- 
teractions in the condensates. For the repulsive bosons 
with point-like interactions K interpolates between 1 in 
the impenetrable Tonks-Girardeu regime and 00 in the 

It is clear that shot noise 

m 



noninteracting regime 18] 



given by the second term in the brackets of Eq. 
is subdominant for large systems as long as if > 1. 
We emphasize that even if shot noise is negligible still 
A4 > A\ as long as K is finite, i.e. there is no true 



long range order in the interfering systems llj. This 
inequality implies that w is finite and the amplitude of 
the interference fringes still fluctuates because of nontriv- 
ial particle-particle correlation functions in each system. 
On the contrary, in systems with short range correlations 
the situation changes and one can not simply ignore shot 
noise. For example, in one-dimensional condensates at fi- 
nite temperature T the particle-particle correlation func- 
tions decay exponentially at long distances [lij ]: 



{a'(x,t)a(x',i)) ~ c 



(16) 



where £t ~ 1/7 1 is the correlation length. Then instead 
of Eq. (TU) we find 



A, 



l + B(pfr) 



l/K-l 



(17) 



with some other nonuniversal coefficients A and B. Note 
that unlike the zero temperature case the second term, 
which corresponds to the shot noise remains finite in the 
thermodynamic limit. At low temperatures p£t 3> 1 its 
contribution remains small but the shot noise becomes 
increasingly important as the temperature grows and £x 
decreases. This result can be readily understood since 
in systems with finite correlation length only particles 
within this length (or the correlation volume in higher 
dimensions) coherently contribute into the interference 
amplitude ll|. The shorter this length the smaller the 
number of such coherent atoms and thus the stronger the 
effects of the shot noise. 

Next we turn to the interference of multiple indepen- 
dent condensates. For simplicity we assume that all the 
condensates are identical and that they form a one di- 
mensional array with a distance d separating the nearest 



neighbors. Then for M condensates the density operator 
n(x, t) assumes the following form in the long time limit: 



n(x, t) « u 2 (x, t) 



M-l 

£« 

j,k=l 



(18) 



where Qj = m(x + dj)/ht and as before we use the fact 
that the Wannier functions corresponding to different 
condensates are identical. We note that using Eq. ([1]) 
as the measure of the interference is not very efficient if 
the number of condensates is large. In particular, one can 
show that the relative width of the distribution of A2 for 
large M is w ~ 1+1/iV, where N is the average number of 
particles per condensate. This width remains finite even 
in the limit of large N. Physically this happens because 
each condensate comes with its own random phase and 
thus there is no constructive interference between them. 
Instead one can define a different interference measure: 



.4: 



q 



dxdx'p(x,t)p(x',t)e iQq(x ~ x,) ~ / dxp{x,t) 



(19) 

where Q q = qQ with q = 1,2, ,,M— 1. We note that in the 
case of multiple condensates there are alternative ways of 
defining such measure. For example in Ref. [l5| Foiling 
et. al. used the quantity C(b) = J dxp(x,t)p(x + b,i). 
The fluctuations of C(b) can be analyzed by methods 
similar to those described in this paper. The evaluation 

of A2 = A is analogous to that of A2 so we only present 
the final answer: 



M(M - 1) 



N 



(20) 



Next we can compute A4, which describes the fluctua- 
tions of A2. The procedure is again quite similar to one 
described above for the case of two condensates. How- 
ever, the actual computations and the resulting expres- 
sion for A4 are quite cumbersome, so we present only the 
two leading terms in the number of the condensates M: 

M 4 N 4 



+ M 3 ( : n 2 : }N 2 - -N 4 + 3N 3 + -N 
6 3 



(21) 

This expression yields the following asymptotic form of 
the relative width of the distribution of the interference 
contrast: 



NVM 



N 2 



5 



(22) 



Note that the width w vanishes with the number of 
the condensates even for the small number of particles 
per condensate. The coefficient multiplying l/y/M in 
Eq. (I22|) weakly depends on the initial state of each con- 
densate. Thus for the interference between bosons in 
number states this coefficient gradually decreases from 
4.2 for N = 1 to 1.6 for N -> 00. 
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A very similar calculation can be done for the case of 
interfering fermions. The only difference between bosons 
and fermions is that all bosonic contributions to A4 come 
with a positive sign and fermionic contributions can bear 
both negative and positive signs. While Eq. (j2"0)) for A 2 
remains intact and holds for fermions as well, the corre- 
sponding expression for A4 gets modified: 



.4 



;(/) „ M*N 4 



Ar 



--N 4 + N 3 + -N 2 
2 3 



(23) 



where the superscript (/) stands for fermions. Cor- 
respondingly the relative width of the distribution be- 
comes: 



,(/) 



NVM 



-N 2 + N + 



(24) 



Note that both Eqs. (|22| and (|24[) have similar scaling 
with the number of the interfering systems. This obser- 
vation leads us to an interesting conclusion that a high 
interference contrast is possible for the case of multiple 
incoherent fermionic sources. This is opposite to the situ- 
ation of two interfering systems, where the fringe contrast 
always remains low. 

We would like to mention that in real experiments de- 
pendence of the interference contrast on M can saturate 
at large M. The reason is that the derivation of Eqs. (j2"2")l 
and (pM|) relied on taking into account all the Fourier com- 
ponents appearing in Eq. (|19p . However, at large wave 
vectors qQ this can become problematic because of fi- 
nite resolution of the apparatus, finite width of Wannier 
functions, not entirely free expansion of the atoms etc. 
Thus realistically the sum over q in Eq. (|19p is always 
bounded by some value q max . Then obviously increas- 
ing the number of interfering systems beyond q ma x will 
not improve the interference contrast. We also point that 
Eq. (TT5|) is valid only in the long time limit provided that 
h 2 t/m > (Md) 2 , or equivalently w > (Md) 2 /wq, where 
w is the width of the condensate after the expansion and 
wo < d is the initial width of each condensate. This con- 
dition is hard to achieve if the number of the interfering 
systems is large. Alternatively one can view the relation 
M* = ^/wwq I d as determining the maximum number of 
the interfering systems beyond which the contrast satu- 
rates. 

In conclusion, we analyzed fluctuations of the ampli- 
tude of interference of independent atomic systems. We 
showed that for the interference of two bosonic systems 
the fluctuations of the fringe amplitude are inversely pro- 
portional to the square root of the number of particles 
in each system. The coefficient of proportionality is non- 



universal; it depends on the details of the wave func- 
tions (or more generally density matrices) of each sys- 
tem. In particular, we found that the fluctuations are 
larger for the case of two interfering condensates with 
broken phases than for the case of two Fock states. We 
generalized our analysis to the interference of multiple 
systems and showed that fluctuations of the fringe con- 
trast vanish as the number of interfering systems (bosonic 
or fermionic) becomes large. 
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